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Abstract 

We describe all countable particle systems on R which have the follow- 
ing three properties: independence, Gaussianity, and stationarity. More 
precisely, we consider particles on the real line starting at the points of 
a Poisson point process with intensity measure m and moving indepen- 
dently of each other according to the law of some Gaussian process ^. We 
describe all pairs (tn,^) generating a stationary particle system, obtaining 
three families of examples. One of these families appeared in connec- 
tion with extremes of independent Gaussian processes in [Z. Kabluchko, 
M. Schlather, L. de Haan, Stationary max-stable fields associated to neg- 
ative definite functions, Ann. Probab. (2009), in press]. 
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1 Introduction 

1.1 Statement of the problem 

We are interested in at most countable systems of particles moving randomly 
on the real line in such a way that the following three requirements are satisfied: 

(Al) The particles are independent of each other. 

(A2) The law describing the motion of each particle is Gaussian and the same 
for all particles. 

(A3) The particles are in an equilibrium. 

The independence stated in requirement (Al) implies that the starting po- 
sitions of particles should be scattered independently over R, which, in more 
rigorous terms, means that they should form a not necessarily homogeneous 
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Poisson point process on K. Requirement (A2) means that the stochastic pro- 
cesses describing the deviations of the particles from their starting positions 
should be Gaussian, having the same law for all particles, and, by require- 
ment (Al), independent of each other. In view of this, the meaning of the first 
two requirements may be described in rigorous terms as follows. 

Let {Ui,i e N} be a Poisson point process on M with intensity measure m. 
We will always assume that m satisfies the following integrability condition: 

e~^^ dm{x) < oo for every £ > 0. (1) 

In most cases of interest, the measure m will be infinite, and so let us agree to 
use N as an index set for the points Ui, even though the case where m is finite 
(and, hence, a.s. only finitely many points Ui exist) is not formally excluded. 

Let € N, be independent copies of a Gaussian process {^{t),t E R"^}. 
We define Vi{t), the position of i-th particle at time t gM.'^ (which we allow to 
be multidimensional), by 

Vi{t) = Ui+^i{t). (2) 

Definition 1. The random collection of functions = {Vi,i G N} will be 
called the independent Gaussian particle system (or simply Gaussian system) 
generated by the pair (m, ^). We use the notation GS{m,^). 

Remark 1. It should be stressed that we do not assume the process ^ to have 
zero mean, which means that we allow for a deterministic component in the 
random motion of particles. In general, it also may happen that ^(0) ^ 0, in 
which case the particles make non-zero jumps immediately after starting at f/j. 

Let us turn to requirement (A3). Given ti,. . . ,tn S R"^, we define a point 
process ^ti,...,t^ on R" by recording the positions of particles at times ti,. . . ,tn- 
That is, we set 

'^t^,...,t„ = {iViiti),...,Viitn)),ieN}. (3) 

The family {^ti,...,t„ : n e N, ti, . . . , f„ e R"^} may be viewed as the family of 
"finite-dimensional distributions" of ^. 

Definition 2. A Gaussian system *p is called stationary if for every n E N, 
every ti, . . . ,tn G K'', and every /i G R'*, we have the following equality of laws 
of point processes on R": 

Vt,+h,...,t„+h = Vt„...,t„- (4) 

1.2 Statement of the main result 

The purpose of this paper is to provide a description of all stationary Gaussian 
systems. Before we can state our main result, we need to recall some definitions. 

A stochastic process {W{t),t G R''} is called stationary if the law of the 
process {W{t + h),t £ R''} does not depend on the choice of /i G R''. A process 
{W{t),t G R**} is said to have stationary increments if the law of {W{t + 
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h) — W{h),t e M"^} does not depend on /i G R'^. Trivially, every stationary 
process has stationary increments. Another family of examples of processes 
with stationary increments is provided by fractional Brownian motions: A zero- 
mean Gaussian process {WK.{t),t G R} is called fractional Brownian motion 
with index k S (0, 2] if 

Gov{w,{h),w^{t2)) = iiir + 1^21" - 1^1 - hr- 

Convention 1 . All stationary processes and processes with stationary increments 
are always supposed to have zero mean. 

A function / : M'^ -> M is called additive if f(ti + ia) = /(^i) + /(i2) for 
every ti,t2 e M'^. Under minor additional assumptions, say, measurability, an 
additive function must be of the form f{t) — {c,t) for some c € M''. 

We will often encounter measures having a density of the form e~^^ with 
respect to the Lebesgue measure. For A e R, we denote by c\ the measure on 
R defined by 

dx 

In particular, Cq is the Lebesgue measure itself. 
The next theorem is our main result. 

Theorem 1. Let S be the set of all pairs (m,^), where m is a measure satis- 
fying H]) and {^{t),t G R"*} is a Gaussian process, with the property that the 
particle system GS{m,^) is stationary. Then 

S = SiUS2US3, (5) 
where the sets Si,S2,S3 are defined as follows: 

1. The set Si consists of all pairs (m,^), where m is an arbitrary measure on 
R satisfying ([T|), and 

{m.i ^^''} = {W{t) + c,t gM!'} 
for some stationary Gaussian process {W{t)^t £ R'*} and some c € R. 

2. The set S2 consists of all pairs (tn, where 

m = aeo and {S,{t),t CM.'^} = {W{t) + f{t) + c,t <=,W^} 

for some a > 0, c € R, a Gaussian process {W{t),t £ W^} with stationary 
increments, and an additive function f : M.'^ —> R. 

3. The set S3 consists of all pairs (tn, where 

m = aeA and {^{t),t eW^} = {W(t) - Xa'^{t)/2 + c,t eW^} 

for some a > 0, A 7^ 0, c G R, and some Gaussian process {W{t), t E W^} 
with stationary increments and variance cr^(t). 
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Although Theorem [T] classifies all pairs (m,^) generating a stationary Gaus- 
sian system, it does not says how to tell whether two given pairs (m', ^'), (m", 
generate equal Gaussian systems or not. This gap will be filled in Section [31 
Now let us state only the following proposition. 

Proposition 1. The decomposition S = where S* = 5i\(iS2U53), is 

disjoint. Pairs belonging to different sets in this decomposition generate different 
in law Gaussian systems. 

1.3 Discussion 

To summarize, we have three disjoint families 5*, ^2, and S3 of stationary 
Gaussian systems. It will be convenient to refer to their representatives as to 
type SI, type S2, and type 53 systems. 

The stationarity of Gaussian systems of type S^ is a rather trivial fact and is 
due to the stationarity of the driving process ^. Somewhat less trivial, but still 
rather appealing is the fact that Gaussian systems of type 52 are stationary. An 
example of a Gaussian system of type ^2 can be obtained by taking m to be the 
Lebesgue measure on M and ^ to be a fractional Brownian motion with a linear 
drift. Surprisingly, the class of stationary Gaussian systems is not exhausted by 
the two "trivial" families SI and ^2: there is one more, non-trivial, family S3. 
An example of a Gaussian system of type S3 can be obtained by taking 

m = ei and {^{t),t e M} = {W^t) - \t\'',te M}, 

where {W^jt e M} is a fractional Brownian motion with index k e (0,2]. For 
K = I, this Gaussian system appeared in [l[ in connection with extremes of 
independent Ornstein-Uhlenbeck processes. For general k G (0,2], the driving 
process WK{t) — |/:|'* appeared in [I4], also in connection with extremes of Gaus- 
sian processes. The results of [l[ were generalized in 0. In particular, it was 
shown in Theorem 2 of [7] that Gaussian systems of type S3 with an additional 
requirement a = l,A=l,c = are stationary. 

The main concern of this paper will be to prove the "difficult" part of The- 
orem [TJ saying that any pair generating a stationary Gaussian system must 
belong to at least one of the three families Si,S2,S3. This will be done in 
Sectional 

At a first sight, it may look that the family ^2 can be included into the family 
S3 by allowing the parameter A to be in the definition of ^3. However, this is 
not the case: the family 52 has an additional "degree of freedom" represented 
by the additive function /. 

In view of particle systems interpretation of Theorem (TJ of special interest 
are Gaussian systems driven by a process ^ satisfying ^(0) = 0. In the next 
corollary we provide a classification of such systems, excluding for convenience 
the non-interesting case in which ^ is a version of the zero process. 

Corollary 1. Let m be a measure satisfying ([T]), and let {^{t),t G M''} be a 
Gaussian process with ^(0) = 0. Assume that for some to, ^{to) is not a.s. 0. 
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Then the particle system GS{m, is stationary iffxn^ ae\ for some a > and 
XeR, and 

j {Wit) + fit), t e R'^} provided that A - 0, 

\ c(t), eGM} — < 

[{Wit) - Xa^it)/2,t eR'^} provided that X ^ 0, 

for some Gaussian process {W{t),t £ R'^} with stationary increments, variance 
(T^(t), and W{0) — 0, and, eventually, an additive function f : R'^ R. 

A natural class of particle systems that has been much studied in the litera- 
ture but is not the subject of the present paper can be obtained by requiring the 
motion of the individual particles to be Markovian rather than Gaussian. More 
precisely, choose the starting positions of the particles accordingly to the law of 
a Poisson point process with cr-finite intensity measure m on some measurable 
space (ft, A), and let the particles move independently of each other accordingly 
to the law of a Markov process on with transition kernel Pix,dy). In this 
setting, the stationarity question has a rather natural solution: It was shown 
in 0] that such particle system is stationary provided that the measure m is 
P- invariant. Various asp ects of Markovian particle systems were studied in 0, 
p. 404], 0, [13], [1], [Hi, see also the references therein. 

Let us also mention that various particle systems driven by Gaussian pro- 
cesses and belonging to type 52 in the above classification were studied in [§], 
14], [^. These papers focus on the behavior of a fixed "tagged" particle in a 



gas of independent Gaussian particles, and their results and methods are rather 
different from ours. 



1.4 Organization of the paper 

Our main result. Theorem [U will be proved in Section [21 The classification 
of stationary Gaussian systems will be completed in Section |3l where we give 
necessary and sufficient conditions for two pairs to generate equal in law systems. 



2 Proof of the main result 
2.1 Notation 

This section will be devoted to the proof of Theorem [TJ We start by introducing 
the notation. We always assume that m is a measure on R satisfying the inte- 
grability condition ([T]), and that {^it),t G R'^} is a Gaussian process. The law 
of the process ^ is uniquely determined by its mean and covariance for which 
we use the notation 

^lit)^Em, r{tut2) = Coviati),ah))- (6) 

Further, we define the variance and the incremental variance of ^ by 

a^t) = Var ^(t), 7(^1, is) = Var [^(ti) - Cih)]- (7) 
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We will often use the identity 

r{h,t2) = i {a\h) + a^{t2) - 7(^1, ^2)) • (8) 

Given ii, . . . , i„ e R'^, the law of the random vector {£,{ti), . . . ,£,{tn)) is denoted 
by ntj_..._t^. 

Let be the Borel cr-algebra of M". For a set B C M" and a; G M, it 

will be convenient to define 

B — X = B ~ (x, . . . , x). 

So, _B — a; is obtained by shifting the set B "diagonally" , in the direction of the 
vector (1, . . . , 1). 

Define '^ti,...,t„, the finite-dimensional distributions of *P, as in ([3|). The 
transformation theory of Poisson point processes, see Proposition 3.8 in (l3l |. 
tells that ^ti....,t„ is a Poisson point process on R" with intensity measure 
nTti,...,t„ that is defined by 

mt„...,tAB)^ fFmti),...,atn))eB^x]dm{x), B e (9) 

In particular, we will often use that vcit — m*nt for every t E M'*, where * denotes 
the convolution of measures. Note that condition ([1]) ensures that mti,....t„{B) 
is finite for every bounded B e ;B(R"). 

We can restate Definition [2] as follows: A Gaussian system *p is stationary 
if for every n&N, every ^i, . . . , t„, /i e M'', and every B e K(M"), 

mt,,...,tAB) ^mt,+h^...,t^+h{B). (10) 

We denote the one-dimensional Gaussian measure with expectation /ig and 
variance CTq by ?l(^o, ctq)- For future reference, let us recall the following formula 
for the Laplace transform of a Gaussian distribution: 

2 2 

if ~ ^{fio,<7l), then E[e«^] = e^°^+^. (11) 
2.2 Two lemmas 

The next two lemmas are standard. We include their proofs only for complete- 
ness. 

Lemma 1. The process W{t) := ^(<) — fi(t) has stationary increments iff for 
all ti,t2,h e M'', 

l{tut2)^j{ti+h,t2 + h). (12) 

Proof. We prove only sufficiency, since the necessity is evident. So, assume 
that (O holds. Let Wh{t) = Wit + h) - W{h). We have 

Cov{Wh{ti),Wh{t2)) = r{ti + h,t2 + h) + r{h,h) - r{h,ti + h) - r{h,t2 + h). 
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Applying ([5]) and taking into account that 7(/i, h) — 0, we obtain 

CoviWhih), Wh{t2)) = hih + h,t2 + h)- j{h, h + h)- -f{h, t2 + h)) . 
By (IT^ , we obtain 

CoY{WH{h),Wh{t2)) - (7(^1, ^2) - 7(0, h) ~ 7(0, t2)) . 

Hence, the law of the process {Wii{t), t G R'^} is independent of h, which proves 
the lemma. □ 

Lemma 2. Let g : M.'^ — > R &e a function satisfying 

g{t2 + h)- giti + h) = g{t2) - g{ti) (13) 
for all ti,t2,h S M.'^. Then the following holds: 

1. The function f[t) := g{t) — g{Q) is additive. 

2. Either g = const or the set of values of g is dense in R. 

Proof. Inserting t2 '■— si, h :— S2, ti :— into (|13p yields /(si + S2) = 
f{si) + /(S2) and proves the first part of the lemma. To prove the second 
part, assume that g is not constant, which means that there is t with f{t) ^ 0. 
A standard inductive argument using the additivity of / gives f{qt) = qf{t) for 
every rational number q. This implies that the set of values of the function /, 
and hence also the set of values of 5, is dense in R. □ 

2.3 Proof of the easy part of Theorem [1] 

In the next three propositions we prove that Gaussian systems of types Si, 82, S3 
are indeed stationary. 

Proposition 2. Let *p ~ GS{m,£,), where (m,^) G Si. Then *p is stationary. 

Proof. By definition of Si, we have the following equality of laws, valid for all 
ti, . . . ,tn, h G M.'^: 

iaum=i^m + hm=i. 

Let B C R" be any Borel set. By ©, we have 



mt,,...,t„(S) = / F[iati),...,atn))eB-z]dmiz) 

JR 

= [ P[(e(ii +h),..., ^{tn + h))eB- z]dm{z) 
Jr 

= ^ti + h,...,tr, + h{B). 

Hence, Eq. (fTU)l holds and Cp is stationary. □ 
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Proposition 3. Let = GS{m,^), where (m,^) G 52- Then *p is stationary. 
Proof. By definition of 52, we have m = aeo for some a > 0, and 

im - c(ti))r=i - m + h)- ah + m^i (w) 

for all ti,...,t„,h€ R'^. Let S C M" be any Borel set. By dH), 

mti,...,t„(5) (15) 

lB-x{yi, ■ ■ ■ ,yn)dnt,,...,t„{yi, ■ ■ ■,yn)dx 
is-(:E+yi)(o,y2 ~yi,---,yn~ yi)dnt^^...,t„{yi, ■ ■ ■,yn)dx 
1b-z(0,?/2 - yi, ...,?;„ - yi)dnf,,...,f„(?/i, . . .,yn)dz 

mu)-atim=i^B-z]dz. 

Taking h G and replacing ti, . . . ,tn hy ti + h, . . . ,tn + h, we get 

mt,+h,...^t„+h{B) - a / P[(e(t, + h)- + /^))r=i e B - (16) 



Comparing P3|) and and taking into account we obtain that ([TI 

holds. Hence, Cp is stationary. □ 

Proposition 4. Let *p = G5(m, ^), where (m,^) G iSs. T/ien *p is stationary. 

Proof. In the particular case a = 1, A = 1, and c — 0, the stationarity of *P 
was proved in Theorem 2 of The general case follows by a straightforward 
application of affine transformations. □ 



2.4 Proof of Theorem [It identifying the driving process ^ 

In Section 12.31 we have shown that 5i U ^2 U ^3 C S. Here we prove the more 
difficult converse inclusion under an additional assumption on the measure m. 
This is stated in the following proposition. 

Proposition 5. Let m be a measure of the form m — at\ +/3eo for some a > 0, 
/3 > 0, A 7^ 0, and let {^(t),t G R"*} be a Gaussian process. Assume that 
= GS'(m,^) is stationary. Then (m,^) G 5iU52U53, where 61,62,33 are as 
in Theorem[l[ 

The proof of Proposition[5]will be based on the following lemma which allows 
to do explicit calculations with measures which are obtained by taking mixtures 
of diagonally shifted and exponentially weighted bivariate normal laws. 
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Lemma 3. Letn be the law of a bivariate Gaussian vector {Xi,X2) with EX^ = 
Hi, Yai Xi — af for i — 1,2, and Yai{Xi — X2) = 7. Let [ be a measure on 
defined for some k by 

1{B) = ( e~'^'n{B - z)dz, B e ^(IR^). (17) 

Then there is a measure l'^'^^ concentrated on the line {{xi,X2) : Xi — Q\ such 
that the following representation holds: 

1{B) ^ J e-'^n^'^^B - z)dz, BeBiR^). (18) 

The Laplace transform ofl^'^\ defined as ^p^'^^u) — e^'^^dl^'^^xi, X2), is given 
by 

ip^'^\u) = exp {(k- u)(/ii + \k(tI) +u{fi2 + ^KCTa) + l-u{u - k)-/ \ . (19) 



Remark 2. Eq.([T9]) shows that the measure [''^■' is a multiple of a Gaussian 
measure with explicitly computable parameters. However, we will not need 
this. 

Remark 3. If the Gaussian measure n has a density, then it is possible to com- 
pute the density of I directly from its definition, Eq. pT]) . However, since n (and 
also I) may fail to have a density, we use a somewhat more complicated represen- 
tation of [ as an exponentially weighted shift of the essentially one-dimensional 
measure [('*^ given in 

Proof of Lemma\^ Define 

[(-^H-B) e^^Us (0,2:2 -xi)dn(xi, 2:2), B€B{M.^). (20) 

By construction, the measure I*-"' is concentrated on the line {(xi, X2) '■ xi = ()\. 
Using transformations similar to those in 0], we obtain 

KB)= [ [ e-"'nB-z{xi,X2)dn{xi,X2)dz (21) 

JR JR2 

e-<''+'''^ e''''nB-(z+x^){^,X2 - xi)dn{xi,X2)dz 
"e'*^ilB-t«(0,X2 - xi)dn{xi,X2)dw. 



R J«? 

Applying (^0)1 to the right-hand side of the above equation, we obtain ^5)1 . 

Now we compute ijj'^'^^u), the Laplace transform of The Laplace trans- 
form of n is defined as 

^(Mi,M2)= / e"i"^+"^"^dn(a;i,a;2). 
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By a two-dimensional analogue of (fTT|) . ip{ui, U2) is given by 

ip{ui,U2) = exp l^iUi + H2U2 + -^{o'lul + 2ruiU2 + crlul)^ , (22) 
where r = Cov(Xi, X2) = {af + o-| - j)/2. It follows from that 
V'(''Hm)= / e''^ie"(="=^-^i)dn(xi,2:2) =?/'(k-w,u). 



The above equation and ([221) yield (1191) after an elementary calculation. □ 

We will also need the following technical lemma. 
Lemma 4. Fix k ^ 0. Let [ be a measure on admitting a decomposition 

1{B)^ [ e-''''[^^\B - z)dz + [ l^°\B-z)dz, B G ^(M^), (23) 



where l^"-* and l^"-' are measures concentrated on the line {{xi,X2) ■ xi — 0}. 
Then the measures and [(0) are uniquely determined. 

Proof. Fix some set A e B{R). For x > 0, let B,j, be a subset of defined by 
= Uj^£[o,a;](A + y). Then implies that 

{{B,) = (^j\-^'d^ [(«)(A)+a;[(")(A). 

The above is valid for every x > 0, and so, {'^'^\A) and \'^'^\A) are determined 
uniquely. □ 

Proof of Proposition\^ Recall from ^ and ^ that ^(•), i7^(-)) and 7(-, •) are 
the expectation, the variance, and the incremental variance of the process ^. We 
start by deducing three statements about /x, tr^, and 7 under various assumptions 
on a, (3, A. 

Claim 1. Assume that a > 0. Then for all ti,t2 G Hi'', 

Mt2)-M^i) = -^(^'(i2)-a2(ti)). (24) 
Proof. The measure mt = m * rij has a density given by the convolution formula 

^""^^^^ = / (ae-^(^-^)+/3)dni(y)=ae-^^ / e^«dnt(2/) + /?. 



Applying (fTTjl to the first term on right-hand side, we obtain 
d-M^) -xx \ ..u^^ , 1 



da; 



ae--'^ exp \^l{t)X + -a' (t)Y \ + /3. (25) 



By stationarity of *p, we must have mtj = tritj for every t\,t2 G M''- This leads 
to (El). □ 
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Before we can prove our next claim, we need to do some preparations. We 
consider the "two-dimensional distributions" of ^p. Take ti,t2 £ K'' and recall 
that ^ti,t2 = {{Vi{ti),Vt{t2))J € N} is a Poisson point process on M^. By 
its intensity measure tritj^tj is given by 

mt,,t,(B) = J{ae-^'=+P)nt,,tAB-x)dx (26) 

= a e~^^nt^^t2{B — x)dx + [3 I ntt,t2{B — x)dx. 

Applying Lemma [3] twice with k = A, rcsp. k = 0, and n — Ufj^tj, we obtain 
two measures on R^, called and tn^^^^, which are concentrated on the line 

{{xi,X2) : a;i = 0} and have the property that for each Borel set B C M^, 

mt,,t2 {B) = a ( e-^^m\'^]t^ {B - x)dx + (3 ( mf^ {B - x)dx. (27) 

JK. Jr 

Claim 2. Assume that a > 0. Then for aU ti,t2,h £ W^, 

7{ti,t2)=-f{ti+h,t2 + h). (28) 

Proof. By stationarity, mti,t2 = mti+hM+h for aU ti,t2,h e R"^. Applying 
Lemma m to the decomposition ((TT)) . we obtain 

Recall that the measures trij^'f^ and tTi|^^^ were constructed by means of 
Lemma[3]and thus have Laplace transforms given by the right-hand side of (jl 
So, we obtain that the expression (considered as a polynomial in u) 

{\-u) L(ti) + ^a2(ti) ) +u(^i{t2) + ^a^{t2)]+lu{u-X)j{ti,t2) 



does not change if we replace ti, t2 by ti + h,t2 + h. Taking into account that 
by Claim [H 

KU) + ^cr^iU) = fi{U + h) + ^a^iU + h), i = 1, 2, 



we arrive at (|28|) . □ 
Claim 3. Assume that /3 > 0. Then for alHi, t2, ^ € R'^, 

^(^2) - ^iitl) = Ai(t2 + h) - fi{h + h) (29) 

and 

l{h,t2)^l{ti+h,t2 + h). (30) 
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Proof. It follows from the decomposition ([?/)) and Lemma 2] that 




■ti+h,t2+h- 



Using the formula for the Laplace transform of tnj^ and ^ti+ht2+h given 
in , we obtain that the expression (considered as a quadratic polynomial in 




remains unchanged if we replace ti , t2 hy ti + h,t2 + h. This yields (|29p and ([50]) . 

□ 

Now we are ready to complete the proof of Proposition \5\ We distinguish 
three cases. 

Case 1. Assume that a > and /3 > 0. We show that in this case, (m, G Si. 
Combining Claim [T] and Claim [3l we obtain 



Since cr^(i) > 0, it follows from Part 2 of Lemma [2] that CT^(i) is a constant 
function. By Claim [U fi{t) is constant as well. Finally, by ^ and Claim [21 



This implies that the Gaussian process W{t) :— ^{t) — ^{t) is stationary. Hence, 



Case 2. Assume that a — Q and /3 > 0. We show that in this case, (m, ^) G 
First of all, note that in this case, m is a multiple of Cq. By Eg. ([30)1 of Claim[3]and 
Lemma [H the process W{t) := £,{t) — fi{t) has stationary increments. Further, 
the function f(t) = n{t) - /Lt(0) is additive by Eq.(l29l) of Claim [3] and Part 1 of 
Lemma [TJ So, we obtain a decomposition ^{t) = W{t) + J{t) + /i(0) implying 
that (m, ^) e ^2. 

Case 3. Assume that a > and /? = 0. We show that in this case, (m, ^) G S^. 
First, we have m = at\. Second, Claim [2] and Lemma [1] show that the process 
W{t) :— ^{t) — ^{t) has stationary increments. It follows from Claim [1] that 



a\ti+h). 



r{ti +h,t2 + h) = - {(T^{ti +h) + (T^{t2 + h) 

^1 (2^2(0) -7(ti,i2)) 
= r(ii,i2). 



-f{h + h,t2 + h)) 



\a^{t)/2 + n{0) + X<T^{0)/2 



where we have set c = ^(0) + Act^(0)/2. Hence, (m,^) e ^3. 



The proof of Proposition [5] is complete. 



□ 
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2.5 Lemmas on convolution equations 

In this section we collect several lemmas on solutions of convolution equations. 
Their proofs are based on explicit calculations with Laplace transforms and on 
the result of [1]. 

Lemma 5. Let xiq — D^(/io, ctq) be a Gaussian measure on M. Let mi, m2 be two 
measures satisfying ([1]) such that 

mi * no = m2 * riQ. (31) 

Then mi = m2 . 

Proof. We assume that (t§ > 0, since otherwise, the statement of the lemma is 
trivial. The density of the measure m^ * no, i = 1, 2, is given by the convolution 
formula 

' f-'^^^.M (32) 



dx V27rcro 
1 



='0 



e ^ / e'oe ^'a dmi{y). 
^27r(To Jr 

Define new measures m'l and m2 by 

^lM^=e z = l,2. (33) 

Let ipyn>,{x) — /jg e^*'(im^(y), i = 1,2, be the Laplace transforms of m'l and 
m'2. Note that by ([1]), (pm'^{x) and Lpm'^{x) are finite for all x e M. We may 
rewrite (15^ as follows: 



d{mi * no) 1 



dx y/2TTao 



e / e^dm-(y) (34) 

Jr 



1 — " , ^tj' 
e ^"a e "0 ip^ 



27rcro 

By (|3ip . the densities of the measures mi * no and m2 * no must be equal. Taking 
into account (p4)) . this yields 



By the uniqueness of the Laplace transform, m'l — m'2- Recalling yields 

that mi = m2. This proves the lemma. □ 

Lemma 6. Let ni — ^{fii,af) and n2 = 0T(/Li2,a-i) be two Gaussian measures 
on M such that al < (t|. Let mi and m2 be two measures satisfying ^ such 
that 

mi * ni = m2 * n2. (35) 
Then mi = m2 * — A^i,^! — ^i)- 
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Proof. We may rewrite as 

mi * ni = (m2 * 9T(/i2 - Mi: o'2 - <^i)) * ^i- 

The proof is completed by applying Lemma [SJ □ 

Lemma 7. Let m be a measure satisfying and let no — 51(/io,o'Q) be a 
Gaussian measure such that for some a > 0, (3 > 0, X ^ 0, 

m * no = aeA + /3eo. (36) 

Then m = ae^^^''o/'^e~^^^°tx + (3tQ. 

Proof. Define a measure mi = ae~^ <^o/^e-^f^o^^ _j_ Then the density of the 
measure mi * no can be computed by means of the convolution formula: 

where the second equality follows from pT|) . Hence, 

m * no = mi * no- 

By Lemma [5l we have m = mi. The proof is completed. □ 

Lemma 8. Let ni — ^{^i^crf) and n2 = DT(/i2,cr|) fee fwo Gaussian measures 
on R such that 7^ cr| . Let m be a measure satisfying ([T]) such that 

m * ni = m * n2. (37) 

Then m = acx + /Jeo /or some a > 0, f3 > and A 7^ 0. 

Proof. By symmetry, we may assume that af < (t| . Then Lemma [S] implies 
that 

m = m*no, (38) 

where no = 0T(/i2 — /xi, cr| — fii). By Theorem 3' of 0], every solution m of 
can be represented as a mixture of exponentials, that is we may write 

^ = / e-^ydpix), 

ay Je 

where p is a finite measure on the set i? = {A e M : e~'^'^dno(x) — 1}. Now, 
in our case the measure no is Gaussian, and so (fTTj) shows that E consists of at 
most two points. One of them is always 0, and the second is denoted by A (if 
E = {0}, let A 7^ be arbitrary). Taking a = p({A}) and /3 — p({0}), we obtain 
m — at\ + f3tQ. This completes the proof. □ 
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2.6 Proof of Theorem [Tt Identifying the measure m 

In this section we complete the proof of the inclusion S C 5i U ^2 U ^3 . Let 
(m, ^) be a pair generating a stationary Gaussian system *p = GS'(m,^). Our 
goal is to show that 

(m,0 e5iU52U53. (39) 

The idea of the proof will be to show, whenever possible, that the measure m is 
of the form ae\ + Pzq and then to apply Proposition O In all other cases, we 
will prove that (tn, ^) G iSi . 

Assume for a moment that ^(0) = and Var^(<o) > for some to e M''. 
Under this restriction, the proof takes the following particularly simple form. 
By stationarity, we have mo = mt„. Using ^(0) = 0, this can be written as 
m = m * rif . Applying to this convolution equation the result of [§| as in the 
proof of Lemma [51 we conclude that m must be of the form acx + /Szq ■ Hence, 
Proposition [5] is applicable and ([M)) is proved. 

In the rest of the section, we are occupied with the proof of Theorem [T] in 
full generality. We distinguish three cases. 

Case 1. Assume that the function is not constant. So, there are ti,<2 G K'' 
such that 

THh) + ^'(^2). (40) 
By stationarity of ^P, we must have trif ^ = trif^ and hence, 

m * OT(//(ti), (T2(ii)) = m * ^{ii(t2), o^{t2)). 

Then Lemma HI which is applicable in view of pp]) . implies that m — m\ + /3eo 
for some a > 0, /? > 0, A 7^ 0. An application of Proposition [5] shows that (|39p 
holds. 

Case 2. Assume that cr^(i) = ct^ > is constant. Take some ti,t2 G and 
fix some -d G [0, 1]. Consider ^txM^ ^ point process on R defined by 

^t,M = {U^ + rnti) + (1 - mih), i e N} , (41) 

where the Ui and the are as in Section [Ol Recalling from ([2]) that Vi{t) = 
Ui + ^i{t), we may rewrite dH]) as 

Vt,M = {^^dti) + (l-m{t2),^eN}. (42) 

By Proposition 3.8 of 13], ^ti,t2 is a Poisson point process whose intensity 
measure vht-^^t2 is given by the formula 

where 

fl{tl,t2)^^^l{tl) + il-^)^I{t2) (43) 
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and 

a^h,t2) = + (1 - df)<7^ + 2d{l - d)r{hM)- (44) 

The stationarity of the particle system *p together with representation (|42p 
shows that for every ti,t2,h S M"^, the point processes ^ti,t2 and ^ti+h,t2+h 
must have the same law. Hence, tTitj^t^ = mt-i^+h.t2+h and consequently, 

m*fn(/i(ii, is), ^2(^1,^2)) = m*D^(/i(ti + h,t2 + h),a'^{ti +h,t2 + h)). (45) 

The proof will be completed after we have considered two subcases. 

Subcase 2a. Assume that for some ti,t2,h g R"^, 

r{h,t2)^r{ti+h,t2 + h). (46) 

Take = 1/2 in the definition of the point process *|5ti,t2- Then (|44p and (|46p 
imply that 

a^{tl,t2) d^ih + h,t2 + h). 

By Lemma m applied to (I45p . the measure m is of the form aCA + Pcq for some 
> 0, > 0, A ^ 0. An application of Proposition [5] shows that ([M)) holds. 

Subcase 2b. Assume that for aU ti,t2, ft. e M"*, 

r(ti,i2)=r(ti+/i,i2 + /i). (47) 

This implies that the process W{t) := £^{t) — ^{t) is stationary. 

If the function fi is constant, then (m, ^) € iSi. Therefore, let us assume that 
fj, is not constant. We will show that this implies that m is a multiple of the 
Lebesgue measure. Let 

G = {g e R : m * Sg = m} 

be the set of "periods" of m, where Sg is the Dirac measure concentrated at g. 
Clearly, G is an additive subgroup of M. 

By stationarity of *P, we have mtj = rritj for every ti, is G R*^- Equivalently, 

m * ^i^J-iti), a^) = m* 0T(M(t2), ct^). 

By Lemma [HI this implies that 

H{h)-fi{t2) e G Vii,i2 eK''. (48) 

Since fi is assumed to be non-constant, Eq. ((48|) implies that G ^ {0}, which 
means that m has a non-trivial period. Of course, this is not sufficient to con- 
clude that m is a multiple of the Lebesgue measure, and so, let us use the sta- 
tionarity of the two-dimensional distributions of *p. Recalling and taking 
into account (|T7)) . we obtain that for every ti, is, h £ M.'^, 

a^(ti,t2) ^a^ih +h,t2 + h). 
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Applying Lemma [H to (|45p . we obtain 

ia) - + h,t2 + h)eG Vti, tj, h £ R'^. 

Recalling a formula for jj, given in we arrive at 

i9 • inih) - - + h) + /i(i2 + /i)) + (^(^2) - ^iit2 + h)) e G. 

Note that this is valid for every d G [0, 1]. Assume that in the above expression, 
d appears with a non-zero coefficient for some ti, ^2, h. Then G contains a non- 
trivial interval, and so, we must have G — M.. In other words, the measure m is 
translation invariant. Since by ([1]), m is finite on bounded intervals, this implies 
that m is a multiple of the Lebesgue measure. 
So, let us assume that for every ti,t2,h e M'', 

- Ai(i2) = +h)- m(<2 + h). (49) 

Recall also that we assume that /i is non-constant. Hence, by Part 2 of Lemma[H 
the set of values of the function /i is dense in M. By ([^5]) . the group G must be 
dense in E. 

We claim that in fact, G = M. To prove this, we need to show that G is 
closed. First of all, the measure m is atomless, since if it would have an atom, 
then the invariance under G would imply that m has a dense set of atoms of 
equal mass, which would contradict ([T|). Now, let gi,g2, ■ ■ ■ be a sequence in 
G converging to some g e M. We claim that g £ G. Indeed, for every interval 
[a, b] C M, we have 

tn([a - g,b- g]) = lim m([a - gi,b- gi\) = lim m([a, b]) = m{[a, b]), 

i — *oo i — >oo 

where the first equality holds since m is atomless, and the second equality follows 
from gi € G. This proves that g E G. Therefore, the group G, being dense and 
closed, must be equal to M. 

The fact that G = R means that the measure m is translation invariant and 
thus, must be a multiple of the Lebesgue measure. Therefore, we can apply 
Proposition [5] which shows that ([55]) holds. 

The proof of Theorem [T] is complete. 

3 Pairs generating equal in law Gaussian sys- 
tems 

In this section we give an answer to the following question: Given two pairs 
(m',^') and (m",^") in S, determine whether GS'(m',^') has the same law as 
G5(m",^"). A first step in this direction is Proposition [T] which we now prove. 

Proof of Proposition [II We will show that Gaussian systems generated by pairs 
belonging to different sets in the decomposition S = S* U 52^ S3 differ by their 
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one-dimensional distributions. If (m,^) G 52, then m = aeo for some a > 0, 
and consequently, mt = tn * nt = aco for every t e M"*. If (m,^) G S3, then 
m = at\ for some a > and A ^ 0. Hence, in this case, trit = m * rif = ae\ 
for some a > 0. Finally, if (m,^) G S^, then rUf is not a multiple of Ca, A G M. 
Otherwise, Lemma[7]would imply that the same is true for m, which contradicts 
the assumption (m,^) G S^. □ 

In the sequel, we concentrate on pairs belonging to the same set in the 
decomposition. Let us call a pair (m, £,) belonging to ^2 or ^3 canonical if 
^(0) = 0. A classification of such pairs was given in Corollary [TJ 

Proposition 6. For every (m,^) G 52 there is a unique canonical pair (tn, ^) G 
S2 generating the same Gaussian system as (m,^). 

Proof. To show the existence, set m = m and £^{t) ^ £,{t) — £,{0)- Applying p3|) 
two times, we obtain that for every B C S(R"), 

mt,,...,tAB) f rmU) - ah))7=i &B- z]dz 

JR 

= a [ miiU) - i{tim=i & B - z]dz 
Jr 

= mt^,...,tAB), 

where mtj,...,*^ are the finite-dimensional intensities of G'5'(m,^), cf. Hence, 
(m, ^) and (m, ^) generate equal in law Gaussian systems. 

We prove the uniqueness part. Let (m, £) be a canonical pair. Then m = ato 
and £{t) = W{t) + f it), see Theorem [Tl We will show that the triple (a, W, /) 
is uniquely determined by the finite-dimensional distributions of *P = GS{m, 

First, we have mt = m * rif = aco for every t G M'^, and so, a is uniquely 
determined. Let us turn to the two-dimensional distributions of By ([5]), we 
have 

mo.t ^ a no,f(B - z)dz. 
By Lemma [3J there is a representation 

n^o,t = a ttIq^j {B — z)dz 

for some measure m'^] concentrated on the line {{xi,X2) G : a;i = 0} and 
having the Laplace transform exp{f{t)u + 1/27(0, t)u^}. By Lemma HJ this 
shows that the two-dimensional distributions of *p determine f{t) and 7(0, t) 
uniquely. To see that 7(0, t) determines the law of W uniquely, recall that 
W^(0) = and hence, we may write the covariance function of W in the form 

r{ti,t2) = i(7(0,ii)+7(0,i2)-7(0,ti-t2)). 
This completes the proof of the uniqueness part. □ 



18 



Proposition 7. For every (tn, ^) G ^3 there is a unique canonical pair (tn,^) G 
S3 generating the same Gaussian system as (tn,^). 

Proof. All necessary ingredients are contained in 0. Take|(i) = ^{t)-^{0) and 
m — m*Sc, where c is as in Theorem [TJ The fact that (m, ^) and (m, ^) generate 
equal Gaussian systems was essentially shown in Proposition 11 of [7|. The 
uniqueness part follows under the additional assumption A = 1 from Remark 24 
of The general case is analogous. □ 

The next proposition gives a necessary and sufRcient condition on two pairs 
belonging to to generate equal in law Gaussian systems. 

Proposition 8. Let (m', and (m", be two pairs, both belonging to S* and 
generating Gaussian systems *P' and Then 

«P' = <p" (50) 

iff the following holds: There is a Gaussian variable Nq whose distribution on 
R is denoted by no and which is independent o/^',^", such that 

m' = m" * no and {C'{t),t G R"^} = {^'{t) + No,t e R"^}, (51) 

or 



m = m 



' * no and {^'{t), teR'^} = {^"(t) + No,teR^}. (52) 



Proof. Introduce the notation /i', r', /i", r", etc. as in Section [Ol By definition 
of the family S^, the functions /i', cr'^, /i", a"^ are constant. Therefore, we 
write, say, /i' instead of fi'{t). We may rewrite ([8]) as follows: 

7'(ii, i2) = 2(a'2 - r'{ti,t2)) and 7"(ii,t2) = 2(a"2 - r"(ti,i2)). (53) 

We start by proving the "if" part of the proposition. Assume for concreteness 
that dniD holds. Then, by (HD, 



SB)^ I m"{ti), e{tn)) eB- z]dm"{z) 

P[(e'(ti) + A^o, ■ • • , Citn) +No)eB- z]dm"{z) 

_ P[(r(ii), • . • , e{tn)) eB-{z + y)]dm"(z)dno(y). 
For every non-negative function / : K*^ ^ K the following formula holds: 
/ f f{z + y)dm"{z)dno{y)^ f f{x)d{m"*no){x). 

Hence, 

(^) = Ijaeiti), . ■ ■,eitn)) eB- x]d{m" * no){x) 

= f P[{^'{ti),...,e{tr.))eB~x]dm'{x) 
Jr 

= K,....tjB). 
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This proves ((501) . 

Now we prove the "only if" part of the proposition. Assume that ((50)) holds. 
Without loss of generality we assume that ct'^ < a"^. Define 

no = 

and let Nq ~ no be a Gaussian variable independent of ^' and We will show 
that dnil) holds. 

We start by proving the first equality in ((5T|l . It follows from ((50]) that 
= m" for all t e K''. Equivalently, 

m'*0T(M',<T'2) = m"* 

Then, by Lemma (HI tn' = m" * Hq. This proves the first equality in ((?T|) . 

We claim that the second equality in ((5T|) follows from the following state- 
ment: for all ti,t2 e M'', 

l'{h,t2)^j"{h,h). (54) 
To see this, set for a moment £^'{t) = (_'(t) + Nq. Then 

E|'(t) = fi' + (m" - //) = /i" = Ee"(i). 

Elementary transformations using ((53| and ((54|) yield 

Cov(|'(ii),|'(t2)) 

= r'(ii,i2) + - = r"(<i,t2) = Cov(r'(ti),r(i2)). 

From now on, we are proving (|54p . We need to consider two cases. 

Case 1. Assume that m' = a'tx +/3eo for some a' > 0, /3 > 0, A 7^ 0. It follows 
from = m" that 

m'*OT(/i',fT") -m"* V,^'"). 

The left-hand side of the above equation is of the form ae\ + /3eo for some a. 
Hence, using Lemma III we conclude that m" = a"t\ + f3to for some a" > 0. 
Let us consider the two-dimensional distributions of By 

m;^.,,(i3) ^a' [ e-^^n;^.,,(S - z)dz + f3 ( n',^,,^{B - z)dz, B e B{^^). 
Jk. Jvl 

Applying Lemma (31 twice, we get two measures ^t^l^ and ^t^\^ concentrated 
on {{xi,X2) £ : Xi — 0} such that the following decomposition is vahd: 

m;^.,,(i3) - a' / e-^^m[[%{B - z)dz + (5 ( m'l^^iB ~ z)dz, B s B{R'). 
Jr Jr 

Furthermore, ipti t^i'^)^ Laplace transform of mi'^\^, is given by 

V'L*,(«) = e'^'(*-*^)"V^ (55) 
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Similar calculations can be done for m"^ f.^ . By ([50)1 , we must have mj^ f.^ 
m'/_^ ■ By Lemma |4l this implies 



m 



'(0) _ ^"(0) 



Comparing the Laplace transforms, we obtain (j54p . 

Case 2. Assume that the condition of Case 1 is not satisfied. We define a point 
process ^Pt^^t^ as in (gT]) and with i? = 1/2: we set 

%^t. = {c^; + e:(ti)/2 + e:(t2)/2, z e n} , 

where {U-, i G N} and f i € N, are the starting points and the driving processes 
of the Gaussian system Then ^t^^^^ is a Poisson point process on R whose 
intensity measure mj^ is given by the formula 



m 



\ ^" ' 2" '2 
A simple calculation using ([55)1 shows that 



m'*fn U', 770- +7:r'(ii,i2) • (56) 



Similar calculations can be done for the pair (m",i^"). By ([50)1 . we must 
have rfij^ — m"^ Denoting these equal measures for a moment by mf^^tj, we 
obtain 

mt,,t, * ^ (0, ]j'{t,,t2)) = mt,,t2 * ^ (0, ]i"{tuh) 



Now assume that ([51]) does not hold for some ti,<2 G IR'^- Then Lemma [5] 
implies that rfit^.tj is of the form ac\ + jSco for some a > 0, /3 > and A ^ 0. 
Further, Lemma [7| applied to (jSG)) yields that m' is of the form a'c\ + 13' for 
some a' > 0, > and A 7^ 0. In fact, the assumption (m',^') G 5* implies 
that even a' > 0, /?' > 0. Hence, we are in the situation of Case 1, which is a 
contradiction. 

The proof of Proposition [8] is complete. □ 



4 Open questions 

We have considered only particles moving on the one-dimensional real line (al- 
though we allowed for a multidimensional time). An interesting question is 
whether it is possible to obtain an analogue of Theorem [1] for particles moving 
in a multidimensional Euclidian space. 

Another problem is to classify all stationary systems of particles driven by 
independent Gaussian processes and starting at the points of an arbitrary point 
process (rather than a Poisson point process). It seems that to gain informa- 
tion from the stationarity of the one-dimensional distributions of such particle 
systems, the results of pj should be used instead of that of [sj. 
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